Singularities in loop quantum cosmology 
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We show that simple scalar held models can give rise to curvature singularities in the effective 
Friedmann dynamics of Loop Quantum Cosmology (LQC). We find singular solutions for spatially 
hat Friedmann-Robertson- Walker cosmologies with a canonical scalar held and a negative exponen- 
tial potential, or with a phantom scalar held and a positive potential. While LQC avoids big bang 
or big rip type singularities, we find sudden singularities where the Hubble rate is bounded, but 
the Ricci curvature scalar diverges. We conclude that the effective equations of LQC are not in 
themselves sufficient to avoid the occurrence of curvature singularities. 
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. The absence of singularities may be taken to be a pre-requisite for any fundamental theory of nature, and in 
CSJ ' particular we expect a quantum theory of gravity to resolve the curvature singularities found in Einstein's theory of 
general relativity. Loop quantum cosmology (LQC) [l[ has been shown to cure the big bang singularity and replace 
q ■ it with a big bounce for Friedmann- Robertson- Walker (FRW) spacetimes with a masslcss canonical scalar field @]. 
The bounce occurs when the matter energy density reaches a Planckian value and quantum gravity effects behave 
repulsively. The results have been generalized by considering the approximate modifications to the Friedmann equation 
implied by LQC. In particular, the models are non-singular for more general forms of the scalar field potential [3(, the 
inclusion of anisotropies in the Bianchi I model [H , the Schwarzschild black hole interior Q , and big rip models with 
constant equation of state Q . These results provide hope that maybe a general singularity resolution theory exists 
O , for the full theory of loop quantum gravity (LQC), of which LQC is a specialized model. 

In this letter we study both canonical and phantom scalar fields with exponential potentials using the effective 
$H | Friedmann equations of LQC. We find a class of singular cosmologies characterized by a bounded Hubble parameter, 
W)i but diverging time derivative of the Hubble parameter occurring in finite proper time. While the energy density is 

bounded, the pressure diverges, a type of singularity known in general relativity as a sudden singularity 
£C} . The energy density and pressure of the scalar field are, respectively, given by 
> . 

P =±-4> 2 +v, P = ±-j 2 -v, (i) 
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where the plus sign indicates a canonical scalar field and the minus corresponds to a phantom field. The local 
conservation of energy-momentum leads to the Klein-Gordon equation governing the dynamics of the scalar field 



g: $ + 3 2ty±£M =0j (2) 



with V(4>) representing the scalar field potential. We will consider exponential potentials of the form 

V = Voe^ , (3) 

where without loss of generality we take A > 0. Scalar fields with exponentials commonly arise in effective theories 
derived from dimensional reductions of higher dimensional models, and appear, for example, in the four-dimensional 
effective theory in the ekpyrotic scenario Q. We have found qualitatively similar results for simple polynomial 
potentials. 

In general relativity the Friedmann equations for a spatially flat, isotropic universe are given by 

ff 2 = yp, H = ~{p + p), (4) 



where H = a/a is the Hubble rate, a the scale factor, and k = 8itG 
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General rclativistic and spatially flat cosmologies with canonical scalar fields and exponential potentials arc always 
singular. With a positive potential 0, [l(| EH, Vq > 0, the scale factor is monotonic with a big bang singularity where 
H ^ oo as a — > at a finite proper time in the past of an expanding universe, or the time reverse (a big crunch 
singularity where H — ► — oo as a — > at a finite time in the future). With a negative potential [11], Vq < 0, there 
exist turning points where H = but these are always maxima of the scale factor where H < 0. Solutions have 
cither a big bang singularity in the past or a big crunch singularity in the future, or both. A phantom scalar field 
with a positive exponential potential, Vq > 0, can avoid a big bang singularity, but suffers a big rip 11311 singularity 
where H — > oo and a — > oo at a finite proper time in the future (or the time reverse in the past) \l4. Il5j. Note that 
one cannot have a phantom scalar field with negative potential as the total energy density must be non-negative in a 
spatially flat Friedmann cosmology. 

To incorporate the quantum effects due to LQC, we consider effective Friedmann equations with corrections of the 
form 0, [| 
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where p c is a constant parameter measuring the loop quantum effects. The value of p c is on the order of the Planck 
density, k~ 4 . Note that these equations are approximations to the true quantum dynamics. They have been validated 
from more rigorous constructions of the quantum dynamics for simpler models with a massless scalar field in 0]. 
Some additional quantum state dependant corrections are suggested in [l6l [l7j , but they have yet to be compared to 
rigorous quantum constructions. Furthermore, the effective Friedmann equation of [TtJ has been derived assuming the 
LQC effects occur in kinetic dominated regimes, which will not be the case for our solutions. The classical GR limit 
can be recovered by letting p c — > oo, whence it is easy to see that the GR Friedmann equations (j4|) are recovered. The 
loop effects do not modify the form of the scalar field equations and thus the Klein-Gordon equation is unchanged. 

The prediction of a bounce in LQC with canonical scalar fields and positive potentials can be understood from 
the modified Friedmann equation ([5]). When the matter energy density reaches the critical value p c , the Hubble rate 
becomes zero triggering the bounce. Therefore, in the case of the canonical scalar field with positive exponential 
potential (Vq > 0), the past big bang or future big crunch singularity of the GR solutions is replaced by a bounce and 
the cosmology becomes non-singular. 

Numerical results for scalar field cosmologies with the loop quantum equations of motion are shown in figure [1] for 
the canonical scalar field with negative potential and figure [2] for the phantom scalar. For the numerics we have taken 
A = 1, Vq = T'* -4 , Pc — k -4 and set k = 1. Initial conditions are chosen such that the scalar field is zero and the 
universe is at a turn-around with H(0) = and energy density equal to p c . This implies a bouncing (re-collapsing) 
phase for the canonical (phantom) case. The qualitative results do not depend on the choice of initial conditions. 

Both the canonical and phantom scalar fields show similar behavior. The solutions are characterized by repeated 
bouncing and collapsing phases with the period of oscillations decreasing to zero in finite proper time. As the period 
decreases, the expanding and collapsing phases become shorter in duration and the scale factor approaches a constant 
at the point when the evolution terminates. The scalar field rolls down (up) the potential for the canonical (phantom) 
case reaching infinity at some final value of t. The Hubble rate is bounded, as expected from the LQC Friedmann 
equation ([5]), while the time derivative of the Hubble rate diverges. Thus the solutions have a "quiescent" [l8[ or 
"sudden" Q singularity where a and H are both bounded but H diverges. The Ricci curvature scalar, R = 6H+12H 2 , 
diverges at a finite proper time. 1 

We note that the case of a phantom scalar field with positive exponential potential in LQC was previously studied 
in (2~l| . Our numerical results are in qualitative agreement, but we find that the solutions are in fact singular. 

In order to verify that the singular nature of the solutions is not a numerical artifact, we can find an approximate 
analytical solution for the regime where a(t) is approximately constant, which becomes increasingly accurate as the 
singularity is approached. The total density, p, remains bounded as the kinetic energy, </> 2 /2, and potential energy, 
V(<p), separately diverge. Setting p = we can use the definition of the energy density (Q]) and solve for cf>(t) to get 



<Kt) = -— In 
An 



(7) 



We note that although the curvature diverges, these singularities are weak in that they may be geodesically complete 19, 2C 
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FIG. 1: Canonical scalar field with negative exponential potential. Plots of a, 0, H and H as a function of time, obtained from 
numerical simulations for k = 1, Vb = — 1, A = 1 and p c = 1. As initial conditions we set 0(0) = 0, 0(0) = y/2(p c — Vb), and 
H(0) — 0. The dashed curve in the plot of 0(t) is the prediction from the approximate analytical solution 
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FIG. 2: Phantom scalar field with positive exponential potential. Plots of a, (/>, H and H as a function of time, obtained from 
numerical simulations, for k = 1, Vb = 1, A = 1 and p c = 1. As initial conditions we use 0(0) = 0, 0(0) = \/2(Vb — p c ), and 
//(0) = 0. The dashed curve in the plot of 0(t) is the prediction from the approximate analytical solution 0. 
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Here to refers to a reference time near the singularity where the scale factor approaches its constant value and 
4>o = 4>(to). Note that this solution applies to both the canonical and phantom scalar fields approaching the singularity. 
In addition, because the energy density is bounded and the pressure grows, the oscillating solutions occur in a regime 
where the kinetic energy of the field approximately equals the potential energy. The analytical solution diverges at a 
critical value of t given by 

Both <f> and <p diverge at this time, and we see that also the pressure and H diverge, indicating a curvature singularity 
at a finite time, in agreement with the numerical results. A comparison of the numerical solution to the analytic 
approximation shows that the two are in close agreement and the approximation becomes better as the critical time is 
approached. This can be seen in figures Q] and [2] where the red dashed line is the analytic solution JZ]) which matches 
well the numerical solution approaching the singularity at t c . Thus the singularity is a genuine feature of the LQC 
equations and not a numerical instability. 

We have shown that the class of exponential potentials studied here gives rise to a singularity in LQC despite the 
quantum corrections in equation ([5]). We can show that this is a consequence of the potential being unbounded. 
The effective Friedmann equation §5§ requires that p is bounded, and for a canonical field, with non-negative kinetic 
energy, the upper bound on p < p c also requires that V < p c , but the potential is not bounded from below by the 
effective Friedmann equation ([S]) . We can write the evolution equation ^ as 



H = -n 2 {p-V) 



1-2^ 



(9) 



Thus for a canonical field with non-negative kinetic energy, H remains bounded if V((f>) is bounded from below. 
Conversely, H is unbounded only if V((j>) is unbounded from below. Similarly for a phantom field with negative 
kinetic energy one can show that H is unbounded only if V{<p) is unbounded from above. 

If we modify the exponential potential ([3|) to introduce a lower bound, the sudden singularity disappears. In figure 
[3] we show numerical results for a canonical scalar field with a potential of the form 

V{<t>) = V Q e XK(t> + V 2 e^ Krt> , (10) 

where Vo < 0, A > 0, as before, but Vi > and A2 > A. These conditions ensure that the potential is bounded from 
below and grows as a positive exponential for large values of <p. In the numerical run, the parameter values were 
chosen to be Vq = k~ 4 , A = 1, V2 = e~ 6 K~ 4 , and A2 = 2, and we have chosen the same initial conditions as the single 
exponential case. With these parameter values the potential minimum is located at 

Inf-^Ws*- 1 . (11) 



«(A 2 -A) V V 2 X 2 



The numerical run indicates that the solution behaves as in the single exponential case with oscillatory expan- 
sion/collapse phases with decreasing period until the potential minimum is reached, after which the period of oscil- 
lations grows. The solution is non-singular as is shown by the boundedness of both the Hubble rate H and its time 
derivative H. The scalar field passes through the potential minimum and rolls up the potential only to turn around 
at a point where the potential is positive which is roughly given by <fi « 6k _1 in figure [H Similar behavior occurs 
in the phantom case for a double exponential potential with Vo > and V2 < where the potential is now bounded 
from above, and the behavior is also non-singular. 

We have shown that singular solutions exist to the LQC effective Friedmann equations for potentials that are not 
bounded from below (above) for a canonical (phantom) scalar field. It is interesting that sudden singularities, which 
may appear to be somewhat contrived in general relativity 0, appear in quite simple models in LQC. While this 
may indicate that generic singularity resolution is not a feature of loop quantum gravity, it may simply indicate that 
the effective equations of LQC break-down in some specific cases. The scalar field Lagrangians which give rise to 
singular behavior also may themselves be regarded as sufficiently pathological that this need not be considered as a 
significant limitation of LQC. 
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FIG. 3: Canonical scalar field with bounded (double exponential) potential. Plots of a, cj>, H and H as a function of time, 
obtained from our numerical simulations, for k = 1, Vo = — 1, A = 1, V2 = e~ 6 , A2 = 2. As initial conditions we set 4>(0) — 0, 
0(0) = - Vb), and H(0) = 0. 
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